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DEFECTS,  CONSTITUTIVE  BEHAVIOR  AND  CONTINUUM  TOUGHNESS 
CONSIDERATIONS  FOR  WELD  INTEGRITY  ANALYSIS 

INTRODUCTION 

The  prediction  of  welded  structural  integrity  is  a  problem  complicated  by  geometry,  inhomo¬ 
geneous  material  composition  and  nonuniform  deformation  fields.  Even  in  the  best  of  circumstances, 
where  optimum  material  quality  and  fabrication  techniques  are  employed,  integrity  prediction  must 
address  these  issues  which  are  not  present  in  many  other  types  of  structural  problems.  For  these  rea¬ 
sons,  concepts  of  material  and  structural  failure  used  to  predict  welded  structural  integrity  should  be 
formulated  in  terms  of  physical  quantities  and  parameters  which  explicitly  address  and  satisfy  the 
analytical  demands  posed  by  welded  structures. 

The  quantitative  relationship  between  fabrication  quality  and  welded  structural  integrity  has  yet 
to  be  determined  in  a  manner  sufficiently  general  to  be  used  in  the  design  process.  Practical  informa¬ 
tion  can.  of  course,  be  obtained  from  weldment  testing.  If  parameters  used  to  describe  weldment 
failure  are  sensitive  to  specimen  size  and  geometry,  however,  the  performance  of  a  weldment  cannot 
be  used  to  reliably  predict  failure  of  a  welded  structure.  Historically,  this  difficulty  has  been  quite 
common. 

The  approach  taken  here  is  based  on  the  explicit  consideration  of  material  fracture  toughness  at 
the  continuum  scale.  At  the  continuum  scale,  the  energy  density  required  to  produce  material  fracture 
is  relatively  insensitive  to  the  particular  degree  of  multiaxialitv  in  the  stress  and  strain  states.  Thus, 
size  and  geometry  effects  are  avoided  using  an  energy  density  material  continuum  fracture  toughness 
definition.  This  further  implies  that  crack  tip  fracture  toughness  parameters  can,  in  principle,  be 
determined  from  analyses  using  energy  density  concepts  to  predict  test  specimen  behavior  which  is 
relevant  for  a  particular  structural  size  and  geometry. 

In  this  investigation,  the  application  of  energy  density  concepts  is  discussed  and  appropriate 
measures  of  continuum  toughness  are  quantitatively  defined.  Interpretation  of  tension  test  data  for  a 
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set  of  short  arc  welds  is  presented  in  a  manner  consistent  with  continuum  plasticity  and  fracture 
toughness.  Finite  element  analysis  results  for  a  welded  T-section  are  presented.  The  predicted  failure 
initiation  load  and  location  is  shown  to  agree  favorably  with  the  experimentally  observed  failure  load. 

CONTINUUM  MATERIAL  TOUGHNESS  CONCEPTS 


A  continuum  volume  of  material  undergoing  deformation  is  characterized  by  its  multiaxial  stress 
state  a,j  and  multiaxial  strain  state  ttJ .  For  ductile  metals,  which  exhibit  inelastic  deformation,  the 
strain  state  is  a  function  of  the  stress  state  and  stress  history.  The  state  of  stress  and  strain  may  be 
related  by  a  constitutive  formulation  which  defines  strain  increments  de^  from  the  current  stress  state. 


The  strain  energy  per  unit  mass  at  a  given  instant  during  deformation,  will  be 
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where  p  is  the  mass  density.  This  energy  density  incorporates  both  stress  and  strain  into  a  fundamen¬ 
tal  quantity  relevant  to  thermodynamic  descriptions  of  material  deformation  and  damage.  The  energy 
density  is  a  scalar  quantity  which  takes  into  account  all  components  of  the  stress  and  strain  tensors  in 
a  physically  consistent  manner.  Failure  of  the  material,  at  the  continuum  scale,  can  be  associated 
with  the  value  of  the  energy  density  at  which  fracture  occurs.  Thus,  the  material  toughness  is  defined 
as 
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where  w(  is  the  critical  strain  energy  density  value.  The  value  of  w,  ,  for  engineering  purposes,  can 
be  taken  as  a  material  constant.  This  makes  w,  independent  of  the  particular  deformation  state  and  is 
a  unifying  feature  absent  in  failure  criteria  based  on  stress  alone  or  strain  alone. 


For  ductile  metals,  the  material  density  varies  only  slightly,  even  over  large  deformations.  For 


this  reason,  it  is  common  to  define  an  energy  per  unit  volume  density 


with  an  associated  critical  value 


>  ... 
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The  energy  per  unit  mass  is  fundamental,  but  the  energy  per  unit  volume  is  equally  appropriate  for 

constant  volume  deformation  processes. 

For  the  case  of  a  uniaxial  true  stress-true  strain  curve,  corresponding  to  a  one-dimensional  state 
of  deformation,  the  critical  energy  density  corresponds  to  the  area  under  the  uniaxial  stress-strain 
curve,  i.e. 

MV  =  Jo  °d(  (7) 

This  representation  is  desirable  for  use  with  traditional  constitutive  formulations  which  rely  on  uniax¬ 
ial  stress-strain  curves. 

For  a  multiaxial  state  of  stress,  each  of  the  six  stress-strain  pairs,  three  normal  and  three  shear, 
must  be  evaluated  and  summed,  i.e. 
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It  should  be  noted  that  one  or  more  individual  terms  in  the  multiaxial  expression  can  be  negative. 
Their  total,  w,  .  must  be  positive,  however. 

A  complete  failure  prediction  should  include  the  location  of  fracture  initiation  within  the  weld 
system  and  the  value  of  the  applied  load  at  which  failure  occurs.  No  a  priori  selection  of  the  fracture 
initiation  site  should  be  made.  Each  constituent  material  (i.e.  base  metal,  weld  metal  and  heat 
affected  zones)  is  a  candidate  failure  site.  The  continuum  toughness  of  each  constituent  material  can 
be  determined  from  simple  material  tensile  tests,  and  the  results  expressed  in  terms  of  the  critical 
energy  density  value  w,  for  each  constituent  material. 


A  maximum  value  of  the  energy  density  will,  in  general,  exist  for  each  constituent  material  in  a 
structure.  A  relative  energy  density  ratio  may  be  defined  by  dividing  the  energy  density  w,  attained 
during  the  loading  process,  by  the  appropriate  critical  energy  density,  w(. .  When  values  of  this  ratio 


are  less  than  unity  fracture  initiation  does  not  occur,  i.e. 

—  <1  n  =  1 . .V  (9) 

w.. 

L  ‘  J  " 

where  N  is  the  total  number  of  constituent  materials  in  the  structure. 

As  the  load  is  increased,  the  structure  will  continue  to  deform.  When  the  relative  energy  den¬ 
sity  ratio  reaches  unity  at  a  site  in  one  of  the  constituent  materials,  say  nc ,  so  that 

=  1  (10) 
and 


<1  n  =  I . 1 ,/»,.  +  1 . N  (11) 

continuum  fracture  occurs  in  material  n(  .  By  normalizing  the  energy  density  throughout  the  structure 
with  the  critical  energy  density  for  each  material,  the  relative  tendency  for  fracture  to  occur  at  dif¬ 
ferent  points  in  the  structure  can  be  quantified. 

WELDED  T-SECTION  GEOMETRY 

'this  investigation  considers  a  welded  T-section  (Fig.  I)  with  flange  and  web  thicknesses  of 
0.2".  The  flange  width  was  2.0".  The  web  height  was  1.8".  End  plates  were  welded  to  the  T- 
section  to  facilitate  testing  of  the  T-section  in  a  tensile  mode. 

The  weld  line  was  located  midway  between  the  ends  of  the  specimen.  The  weld  itself  was 
unsymmetrical.  in  the  sense  that  the  weld  crowns  on  the  top  and  bottom  of  the  flange  were  offset. 
Radiography  of  the  weld  revealed  a  three-dimensional  lack  of  fusion  defect  across  the  flange.  The 
defect  appeared  to  be  offset  with  respect  to  the  weld,  along  a  line  where  the  weld  thickness  was 
approximately  that  of  the  flange.  Also,  the  defect  appeared  to  span  one-half  of  the  flange  thickness 
and  have  an  out  of  plane  dimension  one-eighth  of  the  flange  thickness. 
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For  the  purposes  of  computational  modeling,  the  weld  crown  was  omitted  from  the  model  and 
the  defect  was  located  at  the  center  of  the  weld  material.  This  brought  sufficient  symmetry  to  the 
problem  so  that  only  one-quarter  of  the  flange  needed  to  be  computationally  modeled.  The  defect 
itself  was  idealized  as  a  “slot-shaped”  tunnel  through  the  flange,  with  rounded  ends  (Fig.  2).  This 
interpretation  of  the  radiographic  data  was  consistent  with  the  location  of  the  defect  in  the  weld 
material  and  the  geometry  of  the  defect  as  known  prior  to  the  destructive  testing  of  the  specimen. 

MATERIAL  CHARACTERIZATION 

Data  from  flat  welded  tensile  specimens  was  available  with  which  to  characterize  material 
behavior.  The  specimens  were  0.2  in.  thick  base  metal  plates  which  had  been  welded  together,  in  a 
similar  manner  as  the  T-section.  They  were  cut  into  flat  tensile  specimens  with  the  weld  line  span¬ 
ning  the  specimen  width. 

The  specimens  tended  to  fracture  in  either  the  base  metal  adjacent  to  the  weld  metal,  or  in  the 
weld  metal  itself.  For  the  purposes  of  analysis,  differentiation  between  base  metal  and  weld  metal 
behaviors  were  desired.  No  explicit  treatment  of  the  heat  affected  zone  was  made  given  the  tensile 
specimen  size  and  geometry.  The  occurrence  of  fractures  in  base  metal  adjacent  to  the  weld  included 
its  effects,  in  an  average  sense,  in  the  base  metal  characterization.  Standard  true  stress-true  strain 
curves  were  developed  from  specimen  load-displacement  data.  The  stress-strain  curves  exhibited  sig¬ 
nificant  scatter  but  no  distinguishing  trend  between  base  metal  failure  and  weld  metal  failure  speci¬ 
mens  (Fig.  3).  The  per  cent  reduction  of  area  (%RA),  however,  while  also  exhibiting  scatter,  were 
clearly  divided  into  two  distinct  groups  of  high  %RA  base  metal  failures  and  low  %RA  weld  metal 
failures.  The  base  metal  failures  fell  into  the  range  of  33.0  -  67.0  %RA.  The  weld  metal  failures 
fell  into  the  range  of  0.5  -  6.3  %RA. 

The  marked  difference  between  base  metal  and  weld  metal  ductility  is  not  apparent  from  uniaxial 
considerations  alone.  Accurate  assessments  of  the  respective  material  toughness  require  an  accurate 
uniaxial  true  stress-true  strain  representation.  Such  a  representation  should  quantitatively  reflect  the 


material  ductility  present  in  three-dimensional  states  of  deformation.  To  accomplish  this,  continuum 
uniaxial  true  stress  true-strain  curves  for  the  base  metal  and  weld  metal  (Fig.  4)  were  generated  using 
a  procedure  previously  applied  to  ductile  materials  [1,2].  This  procedure  is  based  on  results  which 
suggest  that  the  actual  continuum  uniaxial  true  stress-true  strain  curve  will  be  more  nonlinear  than  the 
calculated  (or  extrapolated)  standard  true  stress-true  strain  curve  over  the  full  range  of  material  defor¬ 
mation.  This  procedure  was  applied  to  a  bilinear  approximation  of  the  averaged  base  metal  and  weld 
metal  standard  true  stress-true  strain  curves. 

The  critical  strain  ec  of  the  base  metal  at  fracture  was  calculated  based  on  the  average  base 
metal  %RA  since  the  base  metal  was  essentially  defect  free  and  the  %RA  was  influenced  by  weld 
geometry,  some  of  which  might  enhance  apparent  base  metal  toughness,  and  some  of  which  might 
degrade  apparent  base  metal  toughness. 

The  critical  strain  of  the  weld  metal  was  calculated  based  on  the  maximum  weld  metal  *£RA. 
since  the  presence  of  defects  in  the  weld  metal  would  tend  to  reduce  the  apparent  toughness.  The 
maximum  %RA  served  as  a  lower  bound  of  the  material  toughness. 

COMPUTATIONAL  MODELING 

The  ABAQUS  finite  element  code  [3]  was  used  for  the  analysis.  The  finite  element  grid  is 
shown  in  Fig.  5.  Symmetry  considerations,  as  discussed  above,  required  that  only  one  quarter  of  the 
T-frame  specimen  be  modeled.  A  total  of  518  type  C3D20H  and  C3D20  continuum  elements  were 
used.  These  are  20-noded  elements  with  quadratic  displacement  interpolation.  The  C3D20H  ele¬ 
ments  also  include  an  independently  interpolated  linear  hydrostatic  stress.  The  hydrostatic  stress  is 
coupled  to  the  constitutive  relation  using  a  Lagrange  nultiplier.  The  use  of  such  hybrid  elements  in 
regions  of  larger  deformation  near  the  defect  prevented  physically  unrealistic  displacement  constraints 
from  propagating  through  the  grid.  (Such  constraints  can  lead  to  an  artificially  stiff'  response  o! 
standard  elements  for  incompressible  or  nearly  incompressible  deformations.) 


Weld  material  uniaxial  e  stress-true  strain  behavior  was  assigned  to  elements  in  the  denser 
portions  of  the  finite  element  grid  adjacent  to  the  defect  and  to  the  elements  of  the  web  adjacent  to  the 


centerline.  Base  material  behavior  was  assigned  to  the  remainder  of  the  model. 

All  analyses  were  performed  with  full  geometric  nonlinearity  to  account  for  large  strains  and 
large  rotations.  Rik’s  algorithm,  which  is  load-displacement  controlled,  was  used  in  the  analysis  to 
ensure  numerical  stability  since  the  net  load  produced  by  the  specimen  deformation  could  decrease  for 
large  specimen  elongations.  In  the  absence  of  fracture  initiation,  this  would  correspond  to  the  plastic 
limit  load.  Prediction  of  failure  by  this  mode  should  not  be  precluded. 

An  incremental,  rate  independent  plasticity  theory  was  used  for  the  material  constitutive  model 
[4J.  Total  strains  in  the  multiaxial  strain  state  eu  and  uniaxial  strain  state  e  are  linearly  decomposed 
into  elastic  and  plastic  components,  i.e.. 


and 


The  yield  function  /  takes  the  form 


£u  =  eu  +  e‘J 


=  f  f  . 


f(Ojj)  =  o{(P) 


(12) 
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where  a„  and  a  are  multiaxial  and  uniaxial  stress  states,  respectively.  The  associated  flow  rule 
governs  plastic  strain  increments  by  the  relation 

df 


cl  tfj  =  X 


do., 
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In  the  case  of  purely  elastic  behavior  X  -  0.  For  active  yielding. 


X  =  df 


do.. 


>0 


(16) 


Plastic  strain  increments  also  satisfy  a  dissipation  equivalence  condition 

ocl cr  =  ou  clc\) 


and  a  consistency  condition 
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The  von  Mises  yield  function 


/(*«)  =  \  (SUSij)'/2  (19) 

was  employed.  The  deviatoric  stress  s:j  is  defined  as 

1 

sij  =  °ij  ~  j  °kk  ay  .  (20) 

where  they  hydrostatic  component  of  stress  is  a  akk/  3. 

The  finite  element  model  of  the  T-section  was  loaded  by  displacements  applied  to  the  end  of  the 
specimen.  A  static  analysis  was  performed.  The  ratio  of  energy  density  in  the  model  to  the  local 
material  toughness  was  calculated  to  be  high  at  the  center  line  of  the  T-section.  adjacent  to  the  weld 
defect  in  the  weld  metal.  The  element  integration  point  for  Element  1 /Point  25  had  the  greatest 
energy  density  ratio.  The  energy  density  is  plotted  against  the  applied  load  for  this  point  in  Fig.  0 
The  energy  density  is  initially  very  low.  but  as  the  material  begins  to  deform  in  a  nonlinear  fashion 
the  rate  of  change  of  energy  density  with  respect  to  the  applied  load  increases  rapidly.  The  fracture 
of  the  Element  1 /Point  25  was  predicted  to  occur  when  the  energy  density,  w .  reaches  the  critical 
energy  density,  vv(..  for  the  weld  material.  From  this  analysis,  ve  reached  the  4087  lb-in. /in. 3  tough¬ 
ness  value  for  the  weld  metal  at  a  predicted  failure  load  of  30,500  lbs.  The  load-displacement 
response  of  the  T-section  predicted  shows  significant  nonlinearity,  in  part  due  to  yielding  in  both  the 
base  metal  and  weld  metal  (Fig.  7). 

While  the  weld  metal  did  exhibit  nonlinearity  in  the  stress-strain  curve,  its  toughness  value  is 
relatively  low  as  compared  with  the  base  metal  toughness.  It  was  assumed  that  once  fracture  initia¬ 
tion  occurred,  the  resulting  crack  would  propagate  rapidly  as  the  load  bearing  cross-section  of  the 
weld  metal  was  reduced  and  the  load  redistributed  itself  over  the  as  yet  uncracked  portion  of  the  weld 
metal.  For  this  reason,  the  predicted  fracture  initiation  load  was  assumed  to  coincide  with  the  failure 
load  for  the  T-frame.  Of  course,  if  the  weld  material  had  sufficient  toughness,  stable  crack  growth 
would  follow  crack  initiation  and  precede  failure  of  the  T-frame.  Stable  crack  growth  can  K 
modeled,  in  principle,  by  the  same  concept,  i.e..  considering  the  local  material  toughness  and  the 
energy  density  in  the  crack  tip  region. 
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COMPARISON  WITH  EXPERIMENTAL  RESULTS 


As  mentioned  above,  a  laboratory  test  of  the  T-section  was  performed  subsequent  to  the  compu¬ 
tational  analysis.  Conducted  under  structurally  dynamic  loading  conditions,  the  test  produced  gross 
longitudinal  strain  rates  of  the  order  of  10°  sec-1.  This  global  rate  translated  into  local  strain  rates  of 
1G1  sec-1  in  the  vicinity  of  the  defect.  The  local  strain  rate  fell  below  the  strain  rate  sensitivity 
threshold  of  102  sec-1  for  the  weld  metal  and  base  metal,  making  the  predicted  failure  load  from  the 
static  analysis  valid  for  comparison. 


A  T-section  without  a  defect  was  tested  to  provide  a  reference  load  versus  time  response. 
Results  for  this  test  are  shown  in  Fig.  8a.  The  specimen  did  not  fracture  under  the  applied  load  his¬ 
tory.  The  load  signal  rises  to  a  peak  and  subsequently  oscillates  as  transducer  ringing  ensues. 


The  T-section  with  the  flange  defect  was  tested  in  an  identical  manner.  The  results  are  shown 
in  Fig.  8b.  The  specimen  fractured,  and  the  onset  of  the  event  can  be  seen  as  a  distinct  drop  in  the 
load  during  the  period  from  zero  load  to  initial  maximum  load.  The  load  signal  increases  after  this 
drop  and  oscillation  due  to  transducer  ringing  ensures.  This  signal  is  distinct  from  the  defect  free  T- 
section.  The  fracture  initiation  load,  coinciding  with  the  First  peak  of  the  load  history,  was  deter¬ 
mined  using  the  output  from  a  strain  gage  mounted  near  the  weld  region.  The  strain  gage  output  was 
related  to  load  by  a  static  calibration  with  the  load  transducer.  The  fracture  initiation  load  of  the  T- 
section  was  32,000  pounds  and  compared  favorably  with  the  predicted  fracture  initiation  load  of 
30.5000  pounds  obtained  by  the  analysis. 


SUMMARY 


Material  toughness  at  the  continuum  scale  may  be  appropriately  defined  using  the  critical  strain 
energy  per  unit  mass  at  fracture.  For  constant  volume  processes,  the  critical  strain  energy  per  unit 
volume  is  equally  appropriate.  The  value  of  this  critical  energy  density  can  be  determined  from  a 
simple  material  test.  Three-dimensional  aspects  of  the  deformation  and  the  corresponding  influence 
on  the  energy  density  must  be  addressed.  For  uniaxial  material  test  geometries,  the  standard  true 
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stress-true  strain  curve  developed  from  specimen  load-displacement  response  is  not  an  accurate 


description  of  the  material  constitutive  response  for  large  deformations  leading  to  fracture.  A  more 
accurate  continuum  true  stress-true  strain  curve  can  be  developed  from  the  uniaxial  data  and  the  per 
cent  reduction  of  area. 

The  failure  load  of  a  welded  T-section  was  predicted  using  this  methodology.  Radiographic 
examination  of  the  T-section  characterized  a  three-dimensional  lack  of  fusion  defect  present  in  the 
weld  material  across  the  flange.  A  separate  set  of  flat  tensile  coupons  with  transverse  welds  provider! 
information  on  the  constitutive  behavior  and  continuum  toughness  of  the  base  metal  and  weld  metal 
The  predicted  failure  load,  obtained  from  nonlinear  finite  element  analysis,  compared  favorably  with 
the  failure  load  obtained  from  a  laboratory  test  of  the  T-section  performed  subsequent  to  the  numeri¬ 
cal  prediction. 

The  methodology  discussed  in  this  paper  has  been  the  basis  for  other  studies  on  related  types  ol 
weld  integrity  problems  involving  defect  tolerance  in  large  components  [5],  the  influence  of  defects  on 
the  apparent  toughness  of  weld  metal  tensile  specimens  |6]  and  weld  system  performance  |7). 
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Fig.  8  —  T-frame  data  from  laboratory  tests 


